The Amplitude Mode in the Quantum Phase Model 
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We derive the collective low energy excitations of the quantum phase model of interacting lattice 
bosons within the superfluid state using a dynamical variational approach. We recover the well 
known sound (or Goldstone) mode and derive a gapped (Higgs type) mode that was overlooked 
in previous studies of the quantum phase model. This mode is relevant to ultracold atoms in a 
strong optical lattice potential. We predict the signature of the gapped mode in lattice modulation 
experiments and show how it evolves with increasing interaction strength. 
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Interactions can have a dramatic influence on the prop- 
erties of superfluids at low temperatures. In the most 
extreme case, such as lattice bosons at commensurate 
filling, interactions drive a quantum phase transition to 
an insulating phase But even inside the superfluid 
phase, interactions may greatly impact on basic proper- 
ties such as the excitations, with the roton minimum in 
the spectrum of Helium providing a well known exam- 
ple. In this letter, we investigate the effect of increasing 
interaction strength on the collective modes in a super- 
fluid of lattice bosons; the latter are usually modeled 
by the Bose-Hubbard model incorporating both hopping 
(parameter t) and local interactions (U). The question 
has been brought into focus by experiments with ultra- 
cold atoms in optical lattices that were able to control 
the interaction strength and even drive a transition to 
the insulating state More recent experiments have 
probed the excitation spectrum of the superfluid in this 
strongly correlated regime by measuring the energy ab- 
sorption rate in response to periodic lattice modulations. 
The interpretation of the lattice modulation experiments 
is complicated by the large magnitude of the perturba- 
tion and the presence of a confining potential ^3|] ; never- 
theless, the results are very suggestive of the existence of 
a gapped collective mode in the superfluid phase [1, Q , 
which is at the center of interest in the present work. 

The Bose-Hubbard model generates two effective low- 
energy field theories for the complex order-parameter 
field ^p, see Fig.[l] for weak interaction Un ^ t, the effect 
of the lattice can be absorbed in an effective band mass 
and one arrives at the Galilean-invariant Gross-Pitaevskii 
theory (nonlinear Schrodinger equation; n is the mean 
filling per site) . Quite remarkably, close to the superfluid- 
insulator transition at [/ ^ tn, the combined action of the 
lattice and the interaction leads to a Lorentz-invariant 
critical theory (nonlinear Klein Gordon equation) Q. 
The first-order time derivative in the Gross-Pitaevskii 
theory defines the density p = jV'P and hence any density 
mode is bound to the phase degree of freedom, result- 
ing in the unique and well known sound (or Goldstone) 
mode. This differs from the Lorentz-invariant critical 
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FIG. 1: Schematic overview of models describing the su- 
perfluid phase of lattice Bosons. For small (dimensionless) 
interactions u — U/2tnz < the Bose-Hubbard model is 

well approximated by the continuum Gross-Pitaevskii action 
S'gp(V')- Close to the critical interaction u Uc — 2, the 
dynamics is described by the critical theory S'Mott(V')- For 
u > 1/n^ , the Bose-Hubbard model is equivalent to the quan- 
tum phase model, analyzed in detail in this work. 



theory, where the second-order time derivative spoils the 
relation between the order-parameter modulus and the 
density; as a consequence, this theory admits the possi- 
bility of independent amplitude (Higgs) and phase (Gold- 
stone) modes. Naturally, the question poses itself, how 
the Higgs mode emerges upon increasing the interaction 
U. This question should be addressed within a micro- 
scopic theory that explicitly accounts for the interplay 
between interactions and the lattice, as is done by the 
quantum phase model [3, i, a suitable approximation 
of the Bose-Hubbard model for large filling n. In what 
follows, we discuss in more detail the two effective low- 
energy theories for the Bose-Hubbard model and their 
interrelation with the quantum phase model and then 
concentrate on the derivation of the amplitude mode and 
its experimental observation. 



We start with the Bose-Hubbard model in a coherent 
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state formulation where the action is given by 
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here, /i is the chemical potential and /9 = l/k-oT is the 
inverse temperature. The bosonic fields V'i denote the 
amplitude of particles in a Wannier state at site i. 

The diverging correlation length <^Mott c)c 1/ z{uc — u) 
near the (commensurate) superfluid-Mott-insulator tran- 
sition at u = U/2tnz = 2 allows for a description in terms 
of a continuum critical theory C-Matti'^') = {hdr^'l'^ / J z + 
Jz{u, - «)KL,JVVP - IV'n + C/|V'lV2, with z = 2d the 
coordination number and J — 2tn. The emergent (Higgs) 
mode involves a collective oscillation of the amplitude [ipl 
of the order parameter with a frequency that vanishes at 
the transition. Such oscillation of the order parameter is 
accompanied by a local change in the non-condensed frac- 
tion, leaving the local density unchanged. Furthermore, 
this mode is independent of the usual sound mode; this 
can be understood from the fact that the order parameter 
vanishes towards the transition while the density remains 
constant. 

Clearly the gapped amplitude mode is absent in the 
opposite regime, of very weak interactions m <C In 
this case, the length scale set by the interaction is the 
well known healing length i^gp = a^tjUn ^ a. The fact 
that is much larger than the lattice constant a allows 
for a continuum (Gross-Pitaevskii) description £gp (V') — 
V'*fi5rV' + t^nCGp|VV'P + C^IV'l''/2. The effective Galilean 
invariance ensures that at zero temperature there is no 
amplitude mode independent of first sound. 

Comparing the regions of validity of the above coarse 
grained theories, see Fig.[l] we observe that these regimes 
are parametrically separated from each other for n'^ \. 
The absence of a diverging length scale in the intermedi- 
ate regime 1 jr? < u < 1 renders the effect of the lattice 
relevant. As the Galilean invariance is explicitly broken, 
the existence of an amplitude mode cannot be ruled out. 
For large site occupancy, the Hubbard model in the in- 
termediate (and strong) interaction regimes is equivalent 
to the (simpler) quantum phase model 0, Q 
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The conjugate 



zMy, describe the local 



with the Josephson coupling J = 2tn 
operators (pi and 6hj , \(pi , in^ 
phase and deviation from mean filling, respectively. The 
derivation of Hqpm from ([1]) involves an integration over 
density fluctuations under the assumption {5n)/n <^ 1 
[H, which is valid for u ~ and subsequent (re-) 

quantization. Here we restrict our considerations to the 
case of integer filling n e N. 



We analyze the quantum phase model within a dynam- 
ical variational approach |1Q] . which accounts for both 
phase and amplitude degrees of freedom and allows us 
to capture the low-energy physics of a depleted conden- 
sate near the Mott-insulator transition. We first derive 
the static properties in a mean-field approach and then 
include dynamics within a Gaussian approximation. Fi- 
nally, we discuss the response of the system to an exter- 
nal lattice modulation, thereby connecting our findings 
with recent experiments . Before proceeding, we note 
that a gapped excitation closely related to the one con- 
sidered here has been identified by Cazalilla et al. 0] us- 
ing a bosonization approach to a system with a strongly 
anisotropic optical lattice potential. 

Our variational wave function has the Gutzwiller form 
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where jn^) is a particle-number state of the quantum 
phase model and Sni\ni) = {rii — n)\ni). The order- 
parameter in the state given by ^ is ipi = {exp{—iipi)) — 
g-i/4(Tig-i<^i ^ so that (fi and ai determine phase and 
amplitude (fluctuations) of ipi, respectively. The wave- 
function 1^-) has a norm (^'|^') = U^^2nJnJn^ = 
exp[— (mTTCTi)^] and is not properly normalized, 
a consequence of the discreteness of the particle num- 
bers. This is not a problem as long as we stay away from 
the transition so that the particle number fluctuation is 
large (ui ^ 1). We will concentrate on this regime. 
The variational energy Cvar = (^'l-ffqPMl^') is given by 
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and is minimized for tpi — Lp^f = and ai = cr^i = in 
the Mott phase. In the superfluid phase CT,„t ^ 0, leading 
to an order-parameter (or condensate fraction) 
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where W{x) is the Lambert-M^ function [12], cf. 
Fig.lHb). 

If this scheme is carried out blindly all the way to 
strong coupling, the transition to the Mott insulator ap- 
pears as of first (rather than second) order. Besides a lo- 
cal minimum at cr,; =0, which is always present, a second 
minimum first appears at Ui, — 16/e^ ~ 2.16 (spinodal 
point). This spurious first order transition (at ui fa 1.47) 
is due to the failure of the variational wave- function when 
the particle number fluctuations become small (Sfif) 1 
near the transition Q. However, as mentioned earlier we 
will use this approach only away from the critical regime 
of >CMott('0) where ^ should be a good approximation. 
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FIG. 2: (a) Spectra for the phonon (dash-dotted hne) and 
the amphtude mode (dotted Una). Black curves correspond 
to u — U / Jz — 0.25, gray curves to m = 1. (b) Gap Ao- of the 
amphtude mode (dotted hne) and the phonon bandwidth Wf^ 
(dash-dotted hne) vs dimensionless coupUng u — U / Jz in the 
superfluid phase. The gray shaded area corresponds to the 
Mott phase. The solid line denotes the condensate fraction 



Next, we describe the fluctuations above the mean-field 
ground state ^ using a time-dependent variational prin- 
ciple. Following Ref. we aim at an effective action 
(we switch to real time) 



Sell — dt C^if — 



(6) 



While the minimizer of this expression provides the exact 
action in case of unrestricted and properly normalized 
variational states |^'), here, we restrict ourselves to the 
class of variational wavefunctions (cf. [3|) 
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This ansatz generates the Lagrangian 
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Cell = ^ fl(T{E., + hip,<i>, - ^cr, - ^('^ ~ *')^ (8) 

i 

+ e-(i+i6S?-?)/4-.e-(i+i6^?-|)/4-. cos(^, - ipj). 

The new parameters and <i>i in the wave function ([7]) 
assume the role of canonically conjugate fields and allow 
the order parameter's amplitude (ci) and phase (ipi) de- 
grees of freedom to acquire dynamics. The interaction 
terms cx U are linear/quadratic in these fields and the 
coupling between phase (ipi) and amplitude (ci) degrees 
of freedom is only through the hopping term cx J. 

We obtain an effective action in terms of the fields 
ai and (pi by solving the Euler-Lagrange equations for 
Si and $i. The static solution of the Euler-Lagrange 
equations for the ai and ipi fields reproduces the mean- 
field ground state Eq. ^ . Expanding the effective action 
around (cr„f, ipmi) to second order, we find decoupled 
amplitude and phase degrees of freedom in C^ii- We use 
a Legendre transformation [{xi,Xi) (xi,!!^;.), where 



X = a, (fl to quantize the real, classical fields introducing 
ladder operators ([ikj^k'] = (^ick')' 
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Here, = y^h/2ma;UJxik) = B^/uj^k), N is the 
number of sites, and the "masses" are given by mo- = 
fi^\ log |?/'o||/2Jz|V'oP and nicp = fi^/U. We finally obtain 
the quadratic Hamiltonian 
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f^<p{^)0l'Pk + ^f^a{^)'^l^k (10) 



describing phase and amplitude fluctuations above the 
mean-field ground state. The dispersions of the modes 
are given by 



V32| log I Vol P{2 + log IVol [1 + 7(k)]}/i 
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where 7(k) — (2/z) X]i=i cos(k • a;) with lattice vectors 
a; . The amplitude mode is characterized by a finite gap 
Act — hLL!^{0) extending throughout the entire range of 
parameters [l3!], cf. Fig. [21 for small u — > the amplitude 
mode becomes nondispersive, i.e, huja-(k.) w Aa- for all 
k. The phase mode, on the other hand, is gap-less and 
characterized by a sound velocity — ^JlJJz\ij^ a / h. 
This is the Gross-Pitaevskii result, up to the factor l^/ioP 
accounting for the depletion of the condensate at large u. 

Because the sound mode corresponds to a density fluc- 
tuation it can be probed directly by measuring the dy- 
namic structure factor related to the density-density re- 
sponse of the system. This has been done in systems of 
ultracold atoms by using Bragg spectroscopy [1J|. The 
gapped amplitude mode described above is not directly 
accessible to Bragg spectroscopy because it does not in- 
volve a density modulation. Rather, it is excited by 
perturbations that act to modulate the particle number 
variance, or equivalently, the distance to the Mott insu- 
lator phase. Experiments in ultracold atoms have done 
just that By modulating the strength of the optical 
lattice potential, those experiments effectively modulate 
the tunneling which is exponentially sensitive to the lat- 
tice depth. We remark that a similar modulation of the 
coupling was proposed to be relevant in experiments on 
Josephson junction arrays measuring the attenuation of 
ultrasound iJ^]; however, only coupling to the phase de- 
grees of freedom was considered in this context. 

To calculate the response to a lattice modulation, we 
return to the classical theory C^h^ cf. Eq. ([8]), and ex- 
tract the kinetic energy by replacing J ^ J + h, keeping 
track of terms cx h. Going through the re-quantization 
procedure we obtain the kinetic-energy operator 
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FIG. 3: The response function to the lattice modulation 
S^™{u}) within the quadratic theory for different values of 
the interaction strength u in two dimensions. The dark gray 
peak relates to the absorption of a single amplitude mode at 
hcj = Act (for illustrational purposes we give the delta-peak 
a finite width). The area in faint gray corresponds to the 
two-phonon continuum excited at energies below twice the 
phonon bandwidth hui < 2W^ (the logarithmic divergence is 
a density-of-states effect peculiar to two dimensions). Both 
absorption probabilities scale as u at small u. 



T = Ta + which we expand to lowest order in u 



(12) 



The operator To- describes the expected direct coupling of 
the lattice modulations to the amphtude mode at k = 0. 
The other term Tip describes pair-excitations of phase 
modes (phonons) by the lattice modulations. Both oper- 
ators scale as for small interactions. 

Accordingly, the linear response function 5''""(cl') = 
X]„ |("-|7'|0) p(5(?iw ~ hhJno), with \n) the eigenstates of 



the unperturbed Hamiltonian (|T0|) and fujJno the energy 
differences between states |0) and is composed of two 
parts including a single mode peak due to the gapped 
mode and a two-phonon continuum: 



5-"(c.) 



2z^|V.o|4|log|^o|p 



■N5{nw- ^„) + Sf-{Lo). 



Both terms scale as u, see Fig. [3l In two dimen- 
sions, the two-phonon continuum is given by S'^'^{uj) oc 
uuj^K{hujy/l - {huj/2W^)^/W^), with K{x) the com- 
plete elliptic function; going to higher dimensions re- 
quires numerical evaluation. 

Expanding the effective Lagrangian to third order 
in the fields Ci and (pi provides us with the most relevant 
decay channel of the amplitude mode which turns out to 
involve two counter- propagating phonons. We find that 
the small (in dimensions larger than 1) phase space for 
such a process leaves the mode under-damped and damp- 
ing even becomes irrelevant in the limit w — > 0. Finally, 
we comment that the same results can be obtained using 



an RPA type calculation; such an approach allows for a 
systematic improvement of the above results and will be 
discussed in a future publication [l^ . 

In summary, we derived a gapped amplitude (Higgs- 
type) mode showing up in the quantum phase model us- 
ing a dynamical variational approach and discussed its 
relevance in the context of superfluid bosonic atoms ex- 
hibiting strong correlations due to the presence of an op- 
tical lattice. Our analysis demonstrates that this mode 
persists down to weak coupling u « where the 

Gross-Pitaevskii description takes over. We note that 
broken translation invariance (due to the presence of a 
lattice) is crucial for the existence of this mode. This is 
also reflected in the experiment where the coupling to this 
mode is introduced through a modulation of the lattice; 
while experiments ha ve p rovided evidence for the pres- 
ence of such a mode their spectral resolution does 
not yet allow for its detailed analysis. Coupling to this 
mode in a Josephson junction array seems difficult due 
to the rigidity of the coupling parameters; on the other 
hand, the presence of a charge density wave in NbSe2 

MM- 
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While our 



may give access to this mode 
analysis applies to the case of commensurate filling, we 
expect our results to remain valid away from this limit 
as the presence of the amplitude mode is connected with 
a squeezed Gaussian wave function and does not require 
the presence of particle-hole symmetry. 
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